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Abstract. It is shown that there is no standard spiking neural P system that simulates 
Turing machines with less than exponential time and space overheads. The spiking neu- 
ral P systems considered here have a constant number of neurons that is independent 
of the input length. Following this we construct a universal spiking neural P system 
with exhaustive use of rules that simulates Turing machines in linear time and has only 
10 neurons. 



1 Introduction 

Since their inception inside of the last decade P systems [TB] have spawned a variety of hybrid 
systems. One such hybrid, that of spiking neural P systems [3], results from a fusion with 
spiking neural networks. It has been shown that these systems are computationally universal. 
Here the time/space computational complexity of spiking neural P systems is examined. 
We begin by showing that counter machines simulate standard spiking neural P systems with 
linear time and space overheads. Fischer et al. [2] have previously shown that counter machines 
require exponential time and space to simulate Turing machines. Thus it immediately follows 
that there is no spiking neural P system that simulates Turing machines with less than 
exponential time and space overheads. These results are for spiking neural P systems that 
have a constant number of neurons independent of the input length. 

Extended spiking neural P systems with exhaustive use of rules were proved computa- 
tionally universal in [?]. Zhang et al. [18j gave a small universal spiking neural P system with 
exhaustive use of rules (without delay) that has 125 neurons. The technique used to prove 
universality in [1] and [T5] involved simulation of counter machines and thus suffers from an 
exponential time overhead when simulating Turing machines. In an earlier version |10| of the 
work we present here, we gave an extended spiking neural P system with exhaustive use of 
rules that simulates Turing machines in polynomial time and has 18 neurons. Here we improve 
on this result to give an extended spiking neural P system with exhaustive use of rules that 
simulates Turing machines in linear time and has only 10 neurons. 

The brief history of small universal spiking neural P systems is given in Table[TJ Note that, 
to simulate an arbitrary Turing machine that computes in time t, all of the small universal 
spiking neural P systems prior to our results require time that is exponential in t. An arbitrary 
Turing machine that uses space of s is simulated by the universal systems given in [4111118] 
in space that is doubly exponential in s, and by the universal systems given in |3I10I15I1"^ in 
space that is exponential in s. 
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Table 1. Small universal SN P systems. The "simulation time" column gives the overheads 
used by each system we simulating a standard single tape Turing machine, f indicates that 
there is a restriction of the rules as delay is not used and f indicates that a more generalised 
output technique is used. *The 18 neuron system is not explicitly given in [llj; it is however 
mentioned at the end of the paper and is easily derived from the other system presented 
in |11] . Also, its operation and its graph were presented in [12j . 



Chen et al. pT] have shown that with exponential pre-computed resources SAT is solvable 
in constant time with spiking neural P systems. Leporati et al. [7] gave a semi-uniform family 
of extended spiking neural P systems that solve the Subset Sum problem in constant time. 
In later work, Leporati et al. |8] gave a uniform family of maximally parallel spiking neural 
P systems with more general rules that solve the Subset Sum problem in polynomial time. 
All the above solutions to NP-hard problems rely on families of spiking neural P systems. 
Specifically, the size of the problem instance determines the number of neurons in the spiking 
neural P system that solves that particular instance. This is similar to solving problems with 
uniform circuits families where each input size has a specific circuit that solves it. lonescu 
and Drago§ ^5, have shown that spiking neural P systems simulate circuits in linear time. 

In the next two sections we give definitions for spiking neural P systems and counter 
machines and explain the operation of both. Following this, in Section[4l we prove that counter 
machines simulate spiking neural P systems in linear time. Thus proving that there exists no 
universal spiking neural P system that simulates Turing machines in less than exponential 
time. In Section [5] we present our universal spiking neural P system, with exhaustive use of 
rules, that simulates Turing machine in linear time and has only 10 neurons. Finally, we end 
the paper with some discussion and conclusions. 
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2 Spiking neural P systems 



Definition 1 (Spiking neural P systems). A spiking neural P system is a tuple U = 
(O, cTi, (72, • ■ ■ , am, syn, in, out), where: 

1.0 — {s} is the unary alphabet (s is known as a spike), 

2. CTi, (72, • ■ ■ I Cm OLTe neurons, of the form Ui = (n^, 1 ^ i ^ rn, where: 

(a) rii ^ is the initial number of spikes contained in oi, 

(b) Ri is a finite set of rules of the following two forms: 

i. E / s^ ^ s;d, where E is a regular expression over s, b ^ 1 and d ^ 1, 
ii. s*^ — > A;0 where A is the empty word, e ^ 1, and for all E / s^ s;d from Ri 
s'^ ^ L{E) where L{E) is the language defined by E, 

3. syn C {1, 2, • • • , to} x {1, 2, • • • , m} are the set of synapses between neurons, where i ^ j 
for all {i,j) e syn, 

4- in, out G {(7i, (72, • • • , (7m} Ore the input and output neurons respectively. 

In the same manner as in jl5j . spikes are introduced into the system from the environment 
by reading in a binary sequence (or word) w € {0, 1}* via the input neuron ai. The sequence 
w is read from left to right one symbol at each timestep. If the read symbol is 1 then a spike 
enters the input neuron on that timestep. 

A firing rule r = E/s'' — > s; c? is applicable in a neuron a.^ if there are j ^ b spikes in CTj 
and € L{E) where L{E) is the set of words defined by the regular expression E. If, at time 
t, rule r is executed then b spikes are removed from the neuron, and at time t + d — 1 the 
neuron fires. When a neuron at fires a spike is sent to each neuron aj for every synapse (i, j) 
in n. Also, the neuron cr^ remains closed and does not receive spikes until time t + d—1 and 
no other rule may execute in ai until time t + d. We note here that in 2b(i) it is standard 
to have a, d ^ 0. However, we have d ^ 1 as it simplifies explanations throughout the paper. 
This does not effect the operation as the neuron fires at time t + d — 1 instead of t + d. A 
forgeting rule r' = s*^ — > A; is applicable in a neuron ai if there are exactly e spikes in ai. 
If r' is executed then e spikes are removed from the neuron. At each timestep t a rule must 
be applied in each neuron if there is one or more applicable rules at time t. Thus while the 
application of rules in each individual neuron is sequential the neurons operate in parallel 
with each other. 

Note from 2b(i) of Definition [1] that there may be two rules of the form E/s'^ s; d, that 
are applicable in a single neuron at a given time. If this is the case then the next rule to 
execute is chosen non-deterministically. The output is the time between the first and second 
spike in the output neuron a^- 

An extended spiking neural P system [15 has more general rules of the form E/s'^ s^; d, 
where 6 ^ p ^ 0. Note if p = then E/s^ — > s^; d is a forgetting rule. An extended spiking 
neural P system with exhaustive use of rules applies its rules as follows. If a neuron ai 
contains k spikes and the rule E/s'' — > s^; d is applicable, then the neuron ai sends out gp 
spikes after d timesteps leaving u spikes in ai, where k — bg + u, u < b and k,g,u S N. Thus, 
a synapse in a spiking neural P system with exhaustive use of rules may transmit an arbitrary 
number of spikes in a single timestep. In the sequel we allow the input neuron of a system 
with exhaustive use of rules to receive an arbitrary number of spikes in a single timestep. 
This is a generalisation on the input allowed by lonescu et al. [4j. We discuss why we think 
this generalisation is natural for this model at the end of the paper. 

In earlier work [15], Korec's notion of strong universality was adopted for small SN P 
systems. Analogously, some small SN P systems could be described as what Korec refers to 
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as weak universality. However, as we noted in other work [5], it could be considered that 
Korea's notion of strong universality is somewhat arbitrary and we also pointed out some 
inconsistency in his notion of weak universality. Hence, in this work we rely on time/space 
complexity analysis to compare the encodings used by the small SN P system in Table [H 

In the sequel each spike in a spiking neural P system represents a single unit of space. 
The maximum number of spikes in a spiking neural P system at any given timestep during a 
computation is the space used by the system. 

3 Counter machines 

The definition we give for counter machine is similar to that of Fischer et al. [5]. 
Definition 2 (Counter maciiine). 

A counter machine is a tuple C — (z, Cm, Q, qq, q^, /), where z gives the number of counters, 
Cm is the output counter, Q = {qo,qi, ■ ■ ■ ,qh} is the set of states, qa,qh € Q are the initial 
and halt states respectively, E is the input alphabet and f is the transition function 

f : {S xQ X g{i)) {{Y, N} x Q x {INC, DEC, NULL}) 

where g{i) is a binary valued function and ^ i ^ z, Y and N control the movement of 
the input read head, and INC , DEC, and NU LL indicate the operation to carry out on 
counter Ci . 

Each counter Ci stores a natural number value a::. If a; > then g{i) is true and if a; = then 
g{i) is false. The input to the counter machine is read in from an input tape with alphabet 
U. The movement of the scanning head on the input tape is one-way so each input symbol is 
read only once. When a computation begins the scanning head is over the leftmost symbol a 
of the input word aw G S* and the counter machine is in state qo. We give three examples 
below to explain the operation of the transition hmction /. 

— f{a,qj,g(i)) = {Y,qk,INC{h)) move the read head right on the input tape to read the 
next input symbol, change to state qk and increment the value x stored in counter Ci by 
1. 

— f{a,qj,g{i)) — {N, q^, DEC{h)) do not move the read head, change to state qk and 
decrement the value x stored in counter by 1. Note that g{i) must evaluate to true for 
this rule to execute. 

— /(a, qj,g{i)) = {N, qk, NULL) do not move the read head and change to state qk- 

A single application of / is a timestep. Thus in a single timestep only one counter may be 
incremented or decremented by 1. 

Our definition for counter machine, given above, is more restricted than the definition given 
by Fischer [2]. In Fischer's definition INC and DEC may be applied to every counter in the 
machine in a single timestep. Clearly the more general counter machines of Fischer simulate 
our machines with no extra space or time overheads. Fischer has shown that counter machines 
are exponentially slow in terms of computation time as the following theorem illustrates. 

Theorem 1 (Fischer [2j). There is a language L, real-time recognizable by a one-tape TM, 
which is not recognizable by any k-CM in time less than T{n) = 2^^ . 



4 



In Theorem [T] a one-tape TM is an offline Turing macliine with a single read only input 
tape and a single work tape, a k-CM is a counter machine with k counters, n is the input 
length and real-time recognizable means recognizable in n timesteps. For his proof Fischer 
noted that the language L = {waw^ \ w G {0, 1}*}, where is w reversed, is recognisable in 
n timesteps on a one-tape offline Turing machine. He then noted, that time of 2^ is required 
to process input words of length n due to the unary data storage used by the counters of the 
k-CM. Note that Theorem [T] also holds for non-deterministic counter machines as they use 
the same unary storage method. 

4 Non-deterministic counter machines simulate spiking neural P 
systems in linear time 

Theorem 2. Let 11 be a spiking neural P system with m neurons that completes its com- 
putation in time T and space S. Then there is a non- deterministic counter machine Cn 
that simulates the operation of U in time 0{T{xr)^m + Tm?) and space 0{S) where is a 
constant dependant on the rules of 11 . 

Proof idea Before we give the proof of Theorem [2] we give the main idea behind the proof. 
Each neuron <Ti from the spiking neural P system 77 is simulated by a counter a from the 
counter machine Cn- If a neuron ct^ contains y spikes, then the counter will have value y. A 
single synchronous update of all the neurons at a given timestep t is simulated as follows. If 
the number of spikes in a neuron ai is deceasing by b spikes in-order to execute a rule, then 
the value y stored in the simulated neuron Ci is decremented b times using DEC{i) to give 
y — b. This process is repeated for each neuron that executes a rule at time t. If neuron tTj 
fires at time t and has synapses to neurons {a^j , . . . } then for each open neuron ai. in 
{a^j , . . . (Ti^} at time t we increment the simulated neuron Ci- using INC{ij). This process 
is repeated until all firing neurons have been simulated. This simulation of the synchronous 
update of 77 at time t is completed by Cn in constant time. Thus we get the linear time 
bound given in Theorem [5) 

Proof. Let 77 = (O, cti, (T2, • • • , cr„i, syn, m, out) be a spiking neural P system where in — ai 
and out ~ a2- We explain the operation of a non-deterministic counter machine Cn that 
simulates the operation of 77 in time 0{T{xr)^m + Tm^) and space 0{S). 

There are m + 1 counters ci,C2,C3,--- ,Cm,Cm+i in Cn- Each counter Ci emulates the 
activity of a neuron ai. If ai contains y spikes then counter Ci will store the value y. The 
states of the counter machine are used to control which neural rules are simulated in each 
counter and also to synchronise the operations of the simulated neurons (counters). 

Input encoding It is sufficient for Cn to have a binary input tape. The value of the binary 
word w G {1, 0}* that is placed on the terminal to be read into Cn is identical to the binary 
sequence read in from the environment by the input neuron ai. A single symbol is read 
from the terminal at each simulated timestep. The counter ci (the simulated input neuron) 
is incremented only on timesteps when a 1 (a simulated spike) is read. As such at each 
simulated timestep t, a simulated spike is received by ci if and only if a spike is received by 
the input neuron ai . At the start of the computation, before the input is read in, each counter 
simulating ai is incremented Ui times to simulated the Ui spikes in each neuron given by 2(a) 
of Definition [TJ This takes a constant amount of time. 
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Fig. 1. Finite state machine G decides if a particular rule is applicable in a neuron given the 
number of spikes in the neuron at a given time in the computation. Each s represents a spike 
in the neuron. Machine G' keeps track of the movement of spikes into and out of the neuron 
and decides whither or not a particular rule is applicable at each timestep in the computation. 
+s represents a single spike entering the neuron and —s represents a single spike exiting the 
neuron. 

Storing neural rules in the counter machine states Recall from Definition [T] that the 
applicability of a rule in a neuron is dependant on a regular expression over a unary alphabet. 
Let r = i?/s^— >s;dbea rule in neuron cr^. Then there is a finite state machine G that accepts 
language L{E) and thus decides if the number of spikes in ai permits the application of r 
in CTj at a given time in the computation. G is given in Figure [T] If gj is an accept state 
in G then j > b. This ensures that there is enough spikes to execute r. We also place the 
restriction on G that x > b. During a computation we may use G to decide if r is applicable 
in (Ti by passing an s to G each time a spike enters ai. However, G may not give the correct 
result if spikes leave the neuron as it does not record spikes leaving ai. Thus using G we may 
construct a second machine G" such that G' records the movement of spikes going into and 
out of the neuron. G' is construct as follows; G' has all the same states (including accept 
states) and transitions as G along with an extra set of transitions that record spikes leaving 
the neuron. This extra set of transitions are given as follows for each transition on s from a 
state Qi to a state gj in G there is a new transition on — s going from state gi to gj in G' that 
records the removal of a spike from G'. By recording the dynamic movement of spikes, G' is 
able to decide if the number of spikes in ai permits the application of r in ai at each timestep 
during the computation. G' is also given in Figure [TJ Note that forgetting rules ^ A; are 
dependant on simpler regular expressions thus we will not give a machine G' for forgetting 
rules here. 

Let neuron ai have the greatest number I of rules of any neuron in 7J. Thus the applicability 
of rules ri, r2, ■ • ■ , n in ai is decided by the automata G'^, Gj, ■ • ■ ,G[. We record if a rule 
may be simulated in a neuron at any given timestep during the computation by recording 
the current state of its G' automaton (Figure [1]) in the states of the counter machine. There 
are m neuron in 77. Thus each state in our counter machine remembers the current states of 
at most ml different G' automata in order to determine which rules are applicable in each 
neuron at a given time. 

Recall that in each rule of the form r = E/s^ — > s;d that d specifies the number of 
timestep between the removal of b spikes from the neuron and the spiking of the neuron. The 
number of timesteps < d remaining until a neuron will spike is recorded in the states of the 
Cij. Each state in our counter machine remembers at most m different values < d. 
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Algorithm overview Next we explain the operation of C77 by explaining how it simulates 
the synchronous update of all neurons in U at an arbitrary timestep t. The algorithm has 3 
stages. A single iteration of Stage 1 identifies which applicable rule to simulate in a simulated 
open neuron. Then the correct number y of simulated spikes are removed by decrementing the 
counter y times (y = 6 or y = e in 2b of Definition [T]) . Stage 1 is iterated until all simulated 
open neurons have had the correct number of simulated spikes removed. A single iteration of 
Stage 2 identifies all the synapses leaving a firing neuron and increments every counter that 
simulates an open neuron at the end of one of these synapses. Stage 2 is iterated until all firing 
neurons have been simulated by incrementing the appropriate counters. Stage 3 synchronises 
each neuron with the global clock and increments the output counter if necessary. If the entire 
word w has not been read from the input tape the next symbol is read. 

Stage 1. Identify rules to be simulated and remove spikes from neurons Recall that 
d = indicates a neuron is open and the value of d in each neuron is recorded in the states 
of the counter machine. Thus our algorithm begins by determining which rule to simulate in 
counter where ii = min{i \ d = forcri} and the current state of the counter machine 
encodes an accept state for one or more of the G' automata for the rules in ai-^^ at time t. 
If there is more than one rule applicable the counter machine non-deterministically chooses 
which rule to simulate. Let r — E/s'' —> s;dhe the rule that is to be simulated. Using the 
DEC{ii) instruction, counter is decremented b times. With each decrement of the new 
current state of each automaton , G2, • • • , GJ is recorded in the counter machine's current 
state. After b decrements of Ci the simulation of the removal of b spikes from neuron (7^^ is 
complete. Note that the value of d from rule r is recorded in the counter machine state. 

There is a case not covered by the above paragraph. To see this note that in G' in Figured] 
there is a single non-deterministic choice to be made. This choice is at state if a spike is 
being removed (— s). Thus, if one of the automata is in such a state our counter machine 
resolves this be decrementing the counter x times using the DEC instruction. If Ci^ = after 
the counter has been decremented x times then the counter machine simulates state gx-i 
otherwise state gy is simulated. Immediately after this the counter is incremented x — 1 times 
to restore it to the correct value. 

When the simulation of the removal of b spikes from neuron 0-^^ is complete, the above 
process is repeated with counter Ci^ where «2 = minli | ^2 > zi, d = forai} and the current 
state of the counter machine encodes an accept state for one or more of the G' automata for 
the rules in at time t. This process is iterated until every simulated open neuron with an 
applicable rule at time t has had the correct number of simulated spikes removed. 

Stage 2. Simulate spikes This stage of the algorithm begins by simulating spikes traveling 
along synapses of the form j) where ii = min{i \ d = 1 forcri} (if d = 1 the neuron is 
firing). Let {(ji, ji), (ii, J2), • ■ • , {ii,jk)} be the set of synapses leaving ai where j„ < ju+i and 
d ^ I in aj^ at time t (if d ^ 1 the neuron is open and may receive spikes) . Then the following 
sequence of instructions are executed INC(ji), INC(j'2), • ■ • , INC(jfc), thus incrementing any 
counter (simulated neuron) that receives a simulated spike. 

The above process is repeated for synapses of the form {12, j) where 12 = min{i \ 12 > 
ii,d = 1 forai}. This process is iterated until every simulated neuron c,; that is open has 
been incremented once for each spike ai receives at time t. 

Stage 3. Reading input, decrementing d, updating output counter and halting If 

the entire word w has not been read from the input tape then the next symbol is read. If this 
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is the case and the symbol read is a 1 then counter Ci is incremented thus simulating a spike 
being read in by the input neuron. In this stage the state of the counter machine changes 
to record the fact that each k ^ d that records the number of timesteps until a currently 
closed neuron will fire is decremented to /c — 1. If the counter Cm, which simulates the output 
neuron, has spiked only once prior to the simulation of timestep t + 1 then this stage will 
also increment output counter Cm+i- If during the simulation of timestep t counter Cm has 
simulated a spike for the second time in the computation, then the counter machine enters 
the halt state. When the halt state is entered the number stored in counter Cm+i is equal to 
the unary output that is given by time between the first two spikes in cr„i. 

Space analysis The input word on the binary tape of Cn is identical to the length of the 
binary sequence read in by the input neuron of 11. Counters ci to Cm uses the same space as 
neurons ai to am- Counter Cm+i uses the same amount of space as the unary output of the 
computation of 77. Thus Cn simulates U in space of 0{S). 

Time analysis The simulation involves 3 stages. Recall that x > b. Let Xr be the maximum 
value for x of any G' automaton thus Xr is greater than the maximum number of spikes 
deleted in a neuron. 

Stage 1. In order to simulate the deletion of a single spike in the worst case the counter 
will have to be decremented Xr times and incremented Xr — I times as in the special case. 
This is repeated a maximum of 6 < times (where b is the number of spikes removed) . Thus 
a single iteration of Stage 1 take 0{xr^) time. Stage 1 is iterated a maximum of m times per 
simulated timestep giving 0{xr^m) time. 

Stage 2. The maximum number of synapses leaving a neuron i is m. A single spike traveling 
along a neuron is simulated in one step. Stage 2 is iterated a maximum of m times per 
simulated timestep giving O(m^) time. 

Stage 3. Takes a small constant number of steps. 

Thus a single timestep of 77 is simulated by Cn hi 0{{xr)'^m + m?) time and T timesteps 
of 77 are simulated in linear time 0{T{xr)^m + Tni^) by Cn- □ 

The following is an immediate corollary of Theorems [T] and [2] 

Corollary 1. There exist no universal spiking neural P system that simulates Turing ma- 
chines with less than exponential time and space overheads. 

5 A universal spiking neural P system that is both small and time 
efficient 

In this section we construct a universal spiking neural P system that applies exhaustive use 
of rules, has only 10 neurons, and simulates any Turing machine in linear time. 

Theorem 3. Let M be a single tape Turing machine with \A\ symbols and \Q\ states that runs 
in time T. Then there is a universal spiking neural P system 77^/ with exhaustive use of rules 
that simulates the computation of M in time 0{\A\\Q\T) and space 0([2i°S2 r2|QI|A|+2|An]T-) 
and has only 10 neurons. 

If the reader would like to get a quick idea of how our spiking neural P system with 10 
neurons operates they should skip to the algorithm overview in Subsection 15.31 of the proof. 



8 



Proof. We give a spiking neural P system Um that simulates an arbitrary Turing machine 
M in linear time and exponential space. 77m is given by Figure [3] and Tables [2] and [3] The 
algorithm for Um is deterministic and is mainly concerned with the simulation of an arbitrary 
transition rule. Without loss of generality we insist that M always finishes its computation 
with the tape head at the leftmost end of the tape contents. Let M be any single tape Turing 
machine with symbols ai, a2, ■ ■ ■ , a^A\ and states qi, 52, ■ • ■ <1\q\, blank symbol ai, and halt 
state g[Q|. 



5.1 Encoding a configuration of Turing machine M 

Each configuration of M is encoded as three natural numbers using a well known technique. 
A configuration of M is given by the following equation 

Ck — Qr, ■ ■ ■ aiaiai a^^ ■ ■ ■ a_3a_2a-iaoaia2a3 • ■ • flj, aiaiai ■ ■ ■ (1) 

where qr is the current state, each Oi is a tape cell of M and the tape head of M, given by 

an underline, is over ao- Also, tape cells a-x and ay both contain ai, and the cells between 

a-x and ay include all of the cells on M's tape that have either been visited by the tape head 

prior to configuration Ck or contain part of the input to M. 

In the sequel the encoding of object p is given by (p). The tape symbols ai, a2, ■ ■ ■ , ol\a\ 

of M are encoded as (ai) = 1, (a2) = 3, . . . , {a\A\) ~ 2|yl| — 1, respectively, and the states 

gi, (72, ■ • . , 9|Q| are encoded as {qi) = 2 A, ((72) = 4^, . . . , (g|Q|) = 2|(3|A, respectively. The 

contents of each tape cell in configuration Ck is encoded as (aj) = {a) where a is a tape 

symbol of M . The tape contents in Equation ([T]) to the left and right of the tape head are 

X y . 

respectively encoded as the numbers X = ^z'^{a^i) and Y = 'Y^z^{aj) where z — 2"" and 

V = [log2(2|Q| 1^1 + 2|74|)]. Thus the entire configuration Ck is encoded as three natural 
numbers via the equation 

(Ck) - [x, y, {qr) + («,)) (2) 

where (Ck) is the encoding of Ck from Equation ([l} and ai is the symbol being read by the 
tape head in cell aq. 

A transition rule qr,ai,aj, D, qu of M is executed on Ck as follows. If the current state is 
Qr and the tape head is reading the symbol ai in cell ao, aj the write symbol is printed to cell 
ao, the tape head moves one cell to the left to a_i if D = L or one cell to the right to ai if 
D = R, and g„ becomes the new current state. A simulation of transition rule qr,ai, aj, D, g„ 
on the encoded configuration (Ck) from Equation 12]) is given by the equation 

^^^^^^ ^ / (f - (f mod z), zY + (q„> + (f mod z)) 

y[zX + z{aj), ^ - mod z), (g„) + mod z)) 

where configuration Ck+i results from executing a single transition rule on configuration Cfe, 
and (6 mod c) — d where d < c, b = ec + d and &, c,d,e € N. In Equation ^ the top 
case is simulating a left move transition rule and the bottom case is simulating a right move 
transition rule. In the top case, following the left move, the sequence to the right of the tape 
head is longer by 1 tape cell, as cell ao is added to the right sequence. Cell ao is overwritten 
with the write symbol aj and thus we compute zY + z{aj) to simulate cell ao becoming 
part of the right sequence. Also, in the top case the sequence to the left of the tape head is 
getting shorter by 1 tape cell thus we compute — — (— mod z). The rightmost cell of the 
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left sequence a_i is the new tape head location and the tape symbol it contains is encoded as 
mod z). Thus the value mod z) is added to the new encoded current state (qu)- For 
the bottom case, a right move, the sequence to the right gets shorter which is simulated by 
Y — mod z) and the sequence to the left gets longer which is simulated by zX + z{aj). 
The leftmost cell of the right sequence oi is the new tape head location and the tape symbol 
it contains is encoded as mod z). 



5.2 Input to Um 



Here we give an explanation of how the input is read into 11 m ■ We also give a rough outline 
of how the input to Um is encoded in linear time. 

A configuration Ck given by Equation is read into IIm as follows. All the neurons of 
the system initially have no spikes with the exception of ctiq which has 31 spikes. The input 
neuron receives X + 2 spikes at the first timestep ti, Y spikes at time t2, and (g^) + (oi) 
spikes at time ^4. We explain how the system is initialised to encode an initial configuration 
of M by giving the number of spikes in each neuron and the rule that is to be applied in each 
neuron at time t. Thus at time ii we have 

h:a^ = X + 2, s''is^r/s^ s;l, 



where on the left aj = k gives the number k of spikes in neuron aj at time ti and on the 
right is the next rule that is to be applied at time U if there is an applicable rule at that 
time. Thus from Figure [3] when we apply the rule s^(s^)*/s ^ s; 1 in neuron (T5 and the rule 
gSiy'gie _> in neuron aio at time ti we get 

(76, cry, (78, erg = X + 2, s^(s^)*/s A; 0, 

(710 = 15, s^ys^^X;0. 



h:(Ji=X + 2, 
0"?, CTs, (79 = F, 



(s^)7s^A;0, 
s7s4^A;0. 



^4 : (7i = AC, 

(74 = 2, 

cio = 3, 



A;0, 



„3 /„2 



/s s; 1, 



A;0. 
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: ai = X, 

0-4, = (g,.) + (Oi), 
(77,(78,(79 = {qr) + (ai) 



s/s ^ s; log2(2:) + 3. 



Forgetting rules are applied to get rid of superfluous spikes (for example see neurons (77,(73, 
and (7g at time ^2)- Note that (74 is closed at time t2 as there is a delay of 2 on the rule 
(s^(s^)*/s^ s^;2) to be executed in (74. This prevents the Y spikes from entering neuron 
£74 when (75 fires at time ^2- At time ^5 the spiking neural P system has X spikes in (7i, K 
spikes in <T2, and (qr) + {ai) spikes in (74 and erg. Thus at time the spiking neural P system 
encodes an initial configuration of M . 

In this paragraph we will show that given an initial configuration of M it is encoded 
as input to our spiking neural P system in Figure [3] in linear time. In order to do this 
we must compute the three numbers that give (Ck) from Equation [2] in linear time. The 
number X is computed as follows: given a sequence a^xO'-x+i ■ ■ ■ a_2a-i the sequence w = 
(a_^)0'°S2(^)-i(a_^+i)0'°S2(^)-i . . . (a_2)0'°S2(^)-i(a_i)0i°e2(^)-i2 is easily computed in time 
that is linear in x. The spiking neural P system 11 input in Figure [5] takes the sequence w and 
converts it into the X spikes that form part of the input to our system in Figure [3l We give 
a rough idea of how Uinput operates (if the reader wishes to pursue a more detailed view the 
rules for Uinput are to be found in Table |4|). The input neuron of Uinput receives the sequence 
w as a sequence of spikes and no-spikes. On each timestep where (a) is read (a) spikes are 
passed to the input neuron ai , and on each timestep where is read no spikes are passed to the 
input neuron. Thus at timestep ti neuron ai receives {a-x) spikes, and at timestep t2 neurons 
(72, CT3, and CT4 receive {a-x) spikes from ui. Following timestep t2, the number of spikes in 
neurons (72, (73, and (74 double with each timestep. So at timestep <iog2(z)+i the number of 
spikes in each of the neurons (72, (73, and (74 is |(a_^). At timestep iiog2(z)+i neurons 02, 
(73 and (74 also receive {a_x+i) spikes from ai giving a total of z{a-x) + (ci-x+i) spikes in 
each of these neurons at time iiog2(z)+2- Proceeding to time i2iog2(2)+2 neurons a2, (73 and 



£74 have z^{a^x) +'Z(a-x+i) + {a-xj^i) spikes. This process continues until X = ^z*(a_i) is 



computed. The end of the process is signaled when the rightmost number in the sequence is 
read. When this number (2) is read it allows the result to be passed to erg via (75. Following 
this (7g sends X spikes out of the system. Note that prior to this 2 being read only forgetting 
rules are executed in £7g thus preventing any spikes from being sent out of the system. Uinput 
computes X in time a;log2(-2;) + 3. Recall from Section ET] that the value of z is dependant 
on the number of states and symbols in M thus X is computed in time that is linear in a;. 
In a similar manner, the value Y is computed by Uinput in time linear in y. The number 
{lr) + (oj) is computed in constant time. Thus the input {Ck) for 77j\/ is computed in linear 
time. 

5.3 Algorithm overview 

To help simplify the explanation, some of the rules given here differ slightly from those in 
the more detailed simulation that follows this overview. The numbers from Equation ([2|). 
encoding a Turing machine configuration, are stored in the neurons of our system as X, Y 



X 
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input 
i 




output 



Fig. 2. Spiking neural P system Uinput- Each circle is a neuron and each arrow represents 
the direction spikes move along a synapse between a pair of neurons. The rules for Uinput are 
to be found in Tabic [H 

and {qr) + {ui) spikes. Equation ([3]) is implemented in Figure [3] to give a spiking neural P 
system 77m that simulates the transition rules of M. The two values X and Y are stored 
in neurons cti and a2, respectively. If X or y is to be multiplied the spikes that encode X 
or Y are sent down through the network of neurons from either cti or a2 respectively, until 
they reach aiQ. Note in Figure [3] that each neuron from (TTjCTs and erg has incoming synapses 
coming from the other two neurons in ay, as and erg. Thus if <T7,(Ts and ag each contain A'^ 
spikes at time tk, and they each fire sending N spikes, then each of the neurons cry, (Jg and 
CTg will contain 2A'' spikes at time tk+i- Given Y the value zY = T'Y is computed as follows: 
First we calculate 2Y by firing cry, and erg, then 4y by firing a^.a^, and erg again. After v 
timesteps the value zY is computed. zX is computed using the same technique. 

Now, we give the general idea of how the neurons compute ^ — (7- mod z) and 
mod z) from Equation ([3]) (a slightly different strategy is used in the simulation). We begin 
with X spikes in cri. The rule [s'^)* / ^ s; 1 is applied in cri sending ^ spikes to (T4. Following 
this (s^)*s'-^ ^'/s^ — > s^; 1 is applied in 0-4 which sends ^ ~ (7" mod z) to cti leaving 
mod z) spikes in 0-4. The values ^ — mod z) and mod z) are computed in a 
similar manner. 

Finally, using the encoded current state {qr) and the encoded read symbol (a^) the values 
z{aj) and from Equation ([3]) are computed. Using the technique outlined in the first 
paragraph of the algorithm overview the value z((qj.) + (ai)) is computed by sending (q^) + (o^i) 
spikes from erg to crio in Figured Then the rule s^(<«'->+<"'»/s^(<«->+<"'»-<«"> ^ s^<"^>; 1 is 
applied in cio which sends z{aj) spikes out to neurons 174 and erg- This rule uses z({qr) + 
(ai)) — (qu) spikes thus leaving spikes remaining in aiQ. This completes our sketch of how 
IIm in Figure [3] computes the values in Equation ([3|) to simulate a transition rule. A more 
detailed simulation of a transition rule follows. 
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output 



Fig. 3. Universal spiking neural P system Um- Each circle is a neuron and each arrow repre- 
sents the direction spikes move along a synapse between a pair of neurons. The rules for Um 
are to be found in Tables [2] and [3l 



5.4 Simulation of qr, oti, aj, L, (top case of Equation ([3])) 



The simulation of the transition rule begins at time with X spikes in cti, y spikes in (T2, 
ilr) + {<^i) spikes in (T4 and dg, and 1 spike in aiQ. As before we explain the simulation by 
giving the number of spikes in each neuron and the rule that is to be applied in each neuron 
at time t. So at time tk we have 

tk ■ <yi — X, 

^4, = {qr) + (a,), 3^'^^^+^^^^ / s ^ s; 1, 

0-10 = 1, s/s ^ s;log2(z) + 3. 

Thus from Figure [3] when we apply the rule s<'?'-)+("i)/s ^ s; 1 in neurons (T4 and at time 
tk we get 

tk+i : ai = X + {qr) + (a,), ^ log2(^) + 6, 

aio = 1, s/s s;log2(z) + 2. 
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tk+2 ■.ai=X + (qr) + (ai), 
(73 = Y + (qr) + {at), 

if (qr) = ((?|Q|) 
if {qr) ^ {qiQi) 
ar, = Y + (qr) + (at), 
(J6 = Y + (qr) + (ai), 

CTlO = 1, 

tk+3 : CTi = X + {qr) + {ai), 

CT4,CT6 = Y+{qr) + {ai), 
ar, £78, (TQ = Y+ {qr) + {ai), 
CTlO = 1> 



(s")*s<«'->+<"'Vs^s;l, 
s/s ^ s;log2(^) + 1. 

s^(s^)*s<9->+<"'Vs^.s;l, 
s/s s;log2(z). 



In timestep above cr3 the output neuron fires if and only if the encoded current state 
encodes the halt state q\Q\. Recall that when M halts the entire tape contents are to the right 
of the tape head, thus only Y the encoding of the right sequence is sent out of the system. 
Thus the unary output is a number of spikes that encodes the tape contents of M . 

Note that at timestep tk+3 the neuron ar receives Y + {qr) + {ai) spikes from each of the 
two neurons as and ag. Thus at time tk+4 neuron aj contains 2{Y + {qr) + (ai)) spikes. In a 
similar manner erg and ag also receive 2{Y + {qr) + {ai)) spikes at timestep tk+3- The number 
of spikes in each of the neurons ar, as and ag doubles at each timestep between tk+a and 

tfe+log2(^)+2- 

tk+4 : ai = X + {qr) + (a.), s2^(s^)*s<«->+<"'V,s^ ^ .s;log2(^) + 3, 

ar,as,ag = 2{Y + {qr) + {ai)), s^(s^)*s2««'->+<"'»/s ^ ,s; 1, 

0-10 = 1, s/s ^ s;\og2{z) - 1. 

tk+5 : ai = X + {qr) + (a.), s2.(^.),^(g.> + (a,)/^. ^ ^. i^g^(^) ^ 2, 

ar, as, ag = 4(y + {qr) + (a,)), s^(s^)*s^(<3'->+<«*»/s ^ s; 1, 

(710 = 1, s/s ^ s;log2{z) -2. 

tk+6 : C71 = X + {qr) + {ai), s^'is^s^'^^^+^^^^/s^ ^ s; log^C^) + 1, 

(77, £78, £79 = &{Y+{qr) + {a^)), s-(s-)*s8«9'^>+<"-»/s ^ s; 1, 

(710 = 1, s/s^s;log2{z)-3. 



The number of spikes in neurons ar, ag, and £79 continues to double until timestep tfe+iog2(z)+2- 
When neurons £77 and £79 fire at timestep tfe+iog2(z)+2 they send ^{Y + {qr) + {ai)) spikes each 
to neuron crio which has opened at time tk+iog2(z)+2 (for the first time in the transition rule 
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simulation). Thus at time tfe+iog2(2)+3 neuron ctio contains z{Y + (g^.) + (cti)) spikes. 



*fe+log2(z)+2 ■■ (Ti = X + {qr) + {ai), 

2; 

(77,(T8,C79 = -{Y + (qr) + 
flO = 1, 



S^(s^)*5f((9.> + (a.»/s^S.l^ 

s/s s; 1. 



ife+iog2(z)+3 : CTi = X + {qr) + {ai), 
(74, ae = 1, 

<T7, £78, C79 = Z{Y + {qr) + (aj)), 
(Tio = + (gr) + {Oii)), 



s/s — > A; 0, 



Note that {zY mod z^) = and also that z{{qr) + {ai)) < . Thus in neuron ctio at time 
ifc+iog2(^)+3 the rule (s^ y sA{<ir)+{ai)) j _^ . separates the encoding of the right side 
of the tape s^^ and the encoding of the current state and read symbol s^'^'''^+^"*^^. To see 
this note the number of spikes in neurons erg and u\q at time ife+iog2(2)+4- 

The rule ,s^(<9'->+<"'»/.s^«9'->+<"'»-<9">-i s^<"'> ; 1, applied in ctk, at timcstop tfe+iog2(2)+4, 
computes the new encoded current state and the encoded write symbol z{aj). To see 
this note the number of spikes in neurons and ctio at time ifc+iog2(2)+5- Note that neuron 
a\ is preparing to execute the rule s^^(s^)*s^'''^+^"*^/s^ ^ s; 1 at timestep ifc+iog2(z)+6: and 
so at timesteps ffe+iog2(z)+4 and tfe+iog2(^)+5 neuron g\ remains closed. Thus the spikes sent 
out from 0-4 at these times do not enter u\. 

ifc+iog2(.)+4 : = X + {qr) + (a,), s2^(s^)*s<<"->+<«^>/s^ ^ s; 3, 

GA,(J^ = zY, {s'')*/s s;l, 

(TW = Z{{qr) + {ai)), s^((9r> + K»/g^((9.> + {a,»-(g„>-l ^ g^K)-!. 



ife+iog2(z)+5 ■■ ui = X + {qr) + {ai), 

CT2 = ZF, 

aio = (g^) + 1, 



/s^-.s;2, 



(s^)7s^s;l, 



ife+iog2(2)+6 ■■ (Ti = X + {qr) + {ai), 
a2 = zY + z{aj), 
(Tio = {qu) + 1, 



s2^(s^)*s<«'->+<"*Vs^^s;l, 



At time ifc+iog2(z)+7 in neuron 0-4 the rule s^(s^)*s*^^ '""'^ ^^s^ s^: 1 is applied sending 
— — (— mod z) spikes to ai and leaving (— mod z) spikes in a4. At the same time in neuron 
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CTs the rule s^(s^)*s^^ "^""^ ^'^/s'' — > A; is applied leaving only mod z) spikes in 0-5. 



*fe+iog2(z)+7 : 0-1 = (q,.) + (tti), 
(72 = 2:5^ + z{aj), 
X 

(74 = , 

Z 

X 

0-5 = — , 

z 

o-io = (gu) + 1, 



X , 
(74 = — mod z, 

z 

(75 = — mod z, 

z 



Q : <Ti — ( — mod zj, 



(72 ^ zY + z(aj), 

0-4 = (qu) + (— mod z), 

0-6 = {qu) + (— mod z), 

0-7, (78, erg = — mod z, 
z 

c^io = 1, 



„('/u> + (f mod z) 



Jq^) + {§ mod 2:) 



/s s; 1, 
/s -> s; 1, 



s/s ~* s;log2(z) + 3. 



The simulation of the left moving transition rule is now complete. Note that the number of 
spikes in ai, CT2, (74, and at timestep ifc+iog,(2)+9 are the values given by the top case of 
Equation Q and encode the configuration after the left move transition rule. 

The case of when the tape head moves onto a part of the tape that is to the left of a-x+i 
in Equation ([1]) is not covered by the simulation. For example when the tape head is over cell 
a^x+i, then X = z (recall a^^ contains ai). If the tape head moves to the left then from the 
top case of Equation ([3]) the new value for the left sequence is X = 0. Therefore we increase 
the length of X to simulate the infinite blank symbols {ai symbols) to the left as follows. The 
rule s^+(9'-)+("i)/s^ s^;l is applied in ai at time ife+iog2(z)+6- Then at time tfe+iog2(z)+7 
the rule (s^)*/s ^ s; 1 is apphed in 0-4 and the rule s^/s^~^ — > A; is appHed in ur,. Thus at 
time ifc+iog2(z)+8 there are z spikes in ai which simulates another ai symbol to the left, and 
there is 1 spike in (75 to simulate the current read symbol ai. 

We have shown how to simulate an arbitrary left moving transition rule of M. Right 
moving transition rules are also simulated in log2(z)+9 timesteps in a manner similar to that of 
left moving transition rules. Thus a single transition rule of M is simulated by Um in log2(z) + 
9 timesteps. Recall from Section [SH z ~ 2'°S2r2|Q||.4|+2|/i|l ^-j^^g ^^le entire computation of M 
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is simulated in 0{\A\\Q\T) time. From SectionOM is simulated in 0([2i°g2 r2|QI|A|+2|A|l]T) 
space. □ 

While the small universal spiking neural P system in Figure[3]simulates Turing machines with 
a linear time overhead it requires an exponential space overhead. This requirement may be 
shown by proving it is simulated by a counter machine using the same space. However, it is 
not unreasonable to expect efficiency from simple universal systems as many of the simplest 
computationally universal models have polynomial time and space overheads I13I14I17] . 

It was mentioned in Section [2] that we generalised the previous definition of spiking neural 
P systems with exhaustive use of rules to allow the input neuron to receive an arbitrary 
number of spikes in a single timestep. If the synapses of the system can transmit an arbitrary 
number of spikes in a single timestep, then it does not seem unreasonable to allow an arbitrary 
number of spikes to enter the input neuron in a single timestep. If the input is restricted to 
a constant number of spikes, as is the case with earlier spiking neural P systems, then the 
system will remain exponentially slow due to the time required to read the unary input into 
the system. 
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neuron 


rules 




(s^)*g(9r> + (<^.>/g ^ s; 1 if D = R 

s2^(s^)*s<«'-> + <"'Vs^ ^ s;log2(2) + 6 if D=L 
gZ + (,,) + (o,.) _^ s^; 1052(2) + 6 if D=L 
s{9r>+(«.>/g A;0 if D=L 


(72 


(g.).,(9.>+(-.>/5 ^ s; 1 if D=L or {qr) = (giQi) 
s2"(s")*s<«'-> + <°'Vs" s;log2(2) + 6 if D=R 
gZ+(,,) + (o,,) _^ s'-log^iz) + 6 if D=R 




^ A; 0, if (g.)/(g|Q|> 

sVs ^ A;0 


0-4 


(s")7s"»s;l 
sVs^ ^ A;0 

s/s A;0 

S^(S^)*S<T m°dB)/^z ^s^;l 
g(f mod _^ ^.Q 


0-5 


s^(s")7s ^ s;l 

(s")*s<'''-> + <"'>/s s;l 
s"(s")*s(^ "Vs" ^ •^;0 
sVs^"^ ^ A;0 

g(f modz)/g^ g.;^ 


0-6 


s2(s")7s -> A;0 
(s^)7s^s;l 

s/s ^ A;0 

g(£ modz)/_^^^.Q 



Table 2. This table gives the rules in each of the neurons cti to (Tq of IIm- In the rules 
above qr is the current state, ai is the read symbol, aj is the write symbol, D is the move 
direction, and g„ is the next state of some transition rule qr, at, aj, D, q^ of M . Note that 
mod z)), mod z)) e {A} the set of encodings for the symbols of M (see Section lOT) . 
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neuron 


rules 




s^{sY/s A;0 


























s;l 


for all m = 2'', 2 < m ^ z and fe G N 




„{i mod z) /„(— mod z) 


> A;0 




0"lO 


S31/sl6 ^ ^.Q 








A;0 








s7s^ ^ A;0 








sVs^ ^ A;0 








s/s s; log2(2:) + 3 
















^z{(9,) + (<:«,)) /gZ{(q,> + (Q.>)- 


-(9u>- 













Table 3. This table gives the rules in each of the neurons cry to aio of Um- See Table [2] for 
some further explanation. 



neuron 


rules 


cri 


s*/s ^ s; 1 


o"2, 0-3, 0-4 


s*/s ^ s; 1 


0"5 


(s-)*s(-)/,_S;l0g,(^) 




(s")*sVs ^ s;i 


ere 


(s")*s<°>/s ^ A;0 







Table 4. This table gives the rules in each of the neurons of Uinput- 
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